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1 Introduction 



First order phase transitions play an important role in statistical mechanics 
as well as in field theory, see for instance [|T|. For example, the finite tem- 
perature phase transition of Quantum Chromo Dynamics and the symmetry 
restoration phase transition of the weak interaction theory are likely to be 
first order. The value of the interface free energy is a parameter of great 
importance. It determines the non-equilibrium properties of the phase tran- 
sition. In this paper we address the question of how to determine the value of 
the interface free energy using Monte Carlo simulations. For this purpose we 
consider much simpler models namely the 2-d q=10 and q=20 Potts models 
[0], which in statistical mechanics are prototype models of systems with first 
order phase transitions. 

Recently two interesting anal5d;ical results have been obtained for Potts 
models in the first order phase transition region. One is a rigourous theory 
of finite size scaling (FSS) P, ^, which has been proven for large values of 
q, although it is likely to hold for all g > 4. However, it has been found 
numerically that very large lattices are required in order to observe the pre- 
dicted behavior The other result are calculations of the 2-d spin-spin 
correlation length at the infinite volume phase transition point fit [§], 01 • The 
exact value of the order-disorder interface free energy Foa. follows [§]. These 
results could then be compared with previously obtained Monte Carlo (MC) 
results, of which the multicanonical [§, ^ proved to be most accurate. Table 

1 displays the exact interface free energy for q=7, 10 and 20 together with 
MC estimates. Although the overall accuracy of the MC estimates is quite 
satisfactory, it seems that there are uncontrolled systematic errors larger than 
the estimated statistical error. It is thus a challenge to improve the methods 
employed for infinite volume estimates from finite volume numerical calcula- 
tions. These questions are addressed in the next two sections of the paper. 
In section 4 we analyze multicanonical MC data for q=10 and 20. Section 5 
summarizes the final conclusions. 

2 Finite Volume Energy Distribution 

We consider L'^ lattices with periodic boundary conditions. In the 2-d Potts 
model the energy density ranges in steps of in the interval —2 < -E < 0. 
Numerical simulations yield statistical estimates of the canonical probability 
density Pl{P, E). In this paper we are concerned with the shape of Pl{P, E) 
in the vicinity of the transition point in order to extract the interface free 
energy Fo.d.- In the transition region one expects for Pl{I3,E) the existence 
of two distinct maxima p^"-^'" and p™"^''^^ corresponding to the ordered and 
disordered bulk states of the system. These are centred at separated values of 
the energy and Inbetween there will be a value of the energy 
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j^mm ^Y^QTj-Q Pi^{^j3^ E) takes its minimal value p™". States with energies close 
to ^E^*™ correspond to configurations containing two interfaces. 

Following [Q, we can write the partition function for the q-state Potts 
model on a L'^ lattice with periodic boundary conditions in the form 

Zl{P) = e-^'f^hi^^) + qe-^'^m) + 0(e-^^)e-^^(^)^' ; 6 > 0. (1) 

Here are /o(/3) and fd{i^) smooth L independent functions, representing the 
free energy densities of the ordered and disordered phases. Furthermore one 
has = min{/o(/5), At (it both free energies become equal to 

f{Pt)- In the vicinity of Pt they possess the expansions 

pm = PfiPt) + {P- Pt)E^ + ^^~f'^\ ^) + ... (2) 



Pi 



with i = o,d respectively. The Ci and the Ei denote infinite volume specific 
heat and energy densities at Pf In 2-d they have the exact values Eo = 



-1.66425, Ed = -0.96820, - C„ = 0.44763 for q=10 and Eo = -1.82068, 



Ed = -0.62653, Cd-Co = 0.77139 for q=20. The large q expansion of p 
gives Co = 18.06(4) and Co = 5.362(3), in addition numerical simulations 
indicate C„ = 12 - 18 for q=10 and Co = 5.2(2) for q=20 |T3 



We obtain the density of states ^l{E) by inverse Laplace transform. 

^l{E) = — / e'^'^^ZLiO dC. (3) 

We now consider the specific form Zl{(3) = exp{—L'^g{f3)), which is the form 
of the partition function of each of the bulk phases of eq.(p, neglecting terms 
0{e''^^) (we understand that this is not a valid procedure for values of E 
close to E^^'^, see below). It follows 

T ^ /'C+200 

n^i^E) = I e^''^^^-^^^)) dC (4) 



2Txi 



This integral is computed using the steepest descent method. The saddle 
point is at C = Ps{E) which is solution of the equation 



E ^^^^^ 



(5) 



dC 

Ps{E) is the value of P for which E is the internal energy of the system. Since 



= -T'CiPsiE)) (6) 
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is negative the straight path ( = (3s + ix, with x real, is the steepest descent, 
and we obtain the saddle point approximation 

n,{E) ^ , (7) 



The gaussian approximation for Pl{P,E) = yJjb)^ ^^^''Ql{E) is obtained 
expanding E around Ei{P), the internal energy at /3, then 



(8) 

This form corresponds to Binder's ansatz for Pi{l3,E) in the vicinity of a 



first order transition point [|14[. In this approximation the locations of the 
maxima E^""^'^ of Pl{(3,E) show no finite volume dependence. In order to 
find the dependence, one keeps one more term in the expansion of g{l3). Then 
9"k=,siE) ^ g"k=, + {E-Em /'|c=/3//lc=/^), Er"-'' ^ E,^{g"l{2gL% 
This effect is due the asymmetry of the fluctuations inside each bulk phase, 
and is not the one advocated in [p!5| . 

A few remarks on the two Gaussian approximation are appropriate here. 

• The two Gaussian approximation predicts at the transition point (5t a 
ratio Rq = p'^'^^'° j p^"-^^''- = q^JCd/Co- Putting in actual numbers for Co 
and Cd we find -Rio = 10.12 for q=10 and i?2o = 21.4 for q=20. The ratio of 
the weights of the ordered and disordered states is precisely equal to q. This 
follows from eq.(|^) and does not rely on the gaussian approximation. 

• The widths of the two peaks of Pl{i3, E) are predicted, namely 

1 



• Inbetween the two maxima the Gaussian approximation is incorrect. In 
this region mixed phase configurations dominate, this is the subject of the 
next section. 



3 Mixed Phase Configurations 

For energies between the two peaks and large L, the probability to be in a 
pure phase is suppressed like exp[—{E—E)'^L'^/ (2T^Ci)]. Mixed phase config- 
urations, where each phase occupies a macroscopic fraction of the system, be- 
come dominant since they are only suppressed like exp{—ao.d. const. L'^"^), 
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where ao.d. is the order disorder interface tension and const. L*^"^ the interface 
area (interface length for 2-d). Sufficiently close to E^^^ the configurations 
with minimum interface area contains two planar interfaces, closed by peri- 
odic boundary conditions. Closer to the peaks configurations with minimum 
interface are made of a single macroscopic bubble of one phase inside the 
other phase. The shape of the bubble is given by the Wulff construction, it is 
spherical when the correlation length becomes infinite. The interface area of 
a two planar interface configuration is 2L'^~^. In the 2-d case, it dominates 
over the spherical bubble configuration for 

E' < E <E+ (10) 

with 

E- =Eo--{Eo-Ed), (11) 

IT 

E+ = Ed + -{Eo-Ed). (12) 

TT 

Such a behavior has been rigourously shown to occur for the order order 
magnetic transitions of the 2-d Ising model in the broken phase |T^. Nu- 



merically it was observed for the 2-d and 3-d Ising models fl^. Two planar 
interface configurations, with a fraction x of the volume in the ordered phase, 
contribute to the partition function as 

Z^'^'^'ix) oc e-^^'^^°(^) X e-(i--)^'/5/'^(/5) x L^'e"^-''-^'"'. (13) 

The first two factors are contributions from the pure phases. The factor 
LP exp{—Fo,d.L'^~^) comes from interface effects, Fo,d. = 'iao.d. and p = d — 1 



in the capillary wave approximation [|T9 



At the infinite volume transition point Pt the free energy densities fo and 
fd equal f{Pt)- Z'^^^^'^{x) will then asymptotically have the ^-independent 
form 

Z-™-^ |^=^,oc e-^'^*^(^*) X LP e-^°-^-^'"\ (14) 

It means that in the region around E'™" the probability distribution 
PLif^uE) = oc LP exp(-Fo,d,L'^-i) is fiat. This formula could be 

used to determine the interface tension. However, we found it more con- 
venient to use the ratio ^pp^^^ppM/pmrn ^ L-p+'^/^ exp{Fo.dL'^-^), and 
introduce the finite volume estimator 

2 pmax ,o pmax ,d 

^(i) = ^ ''(p.Jy I l^=ft (15) 

- + + + + + ^ + (16) 
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It is convenient to consider PL{Pt,E) because one can check whether it has 
a flat part as it should. If it does not, the value of Fo,d. extracted from eq. 
(|16D is questionable. However this equation, with L'^{P — (3t) dependent a^'s, 
holds asymptotically for all /3 in a oc l/L'^ neighborhood of (3t (e.g. at the 
equal heights point like in 0). 



4 Numerical Results 

The simulations have been performed with the multicanonical algorithm [^j, 
for more information we refer to ||^. For q=10 we have simulated lattices 
in the range from 12^ to 100^, for q=20 from 16^ to 70^. The statistics 
for q=10 is given in |0, for q=20 in and additional q=20 statistics is 
collected in Table 2. Errors quoted in this paper are bias corrected jackknife 
estimates. The simulations have been performed on the Saclay Cray X-MP, 
on RISC stations of the SCRI workstation cluster and on an ALPHA station 
at Bielefeld University. 

ror the determmation of the quantities Fj^ i , F^ and Fj^ 

we adopt the following procedure. After locating approximate values we fit 
the probability function FilPt, E) in the vicinity of the maxima on an energy 
interval corresponding to a decrease of FL{j3t,E) by at most a factor 1/e of 
its maximal value. We use the cubic form 

/ 771 _ rpmax,i\2 

In FM, E) = const + ^ — ^ ^ + a{E - ^r"'')'- (17) 

We have checked that the fits gave acceptable x^/d-o.f values. Figure |I| 
displays an example for the q=20 model on a 50^ lattice. It can be seen 
that the probability distribution in vicinity of the maximum shows a sizable 
asymmetry, demonstrating that the cubic term is needed in order to give a 
good description of the data. For the determination of the quantities P|"*" 
and E^^"- we found it sufficient to use a quadratic fit analogue to eq.(|T7|) 
with a = 0. An example of such a fit is displayed in Figure |. 

Using eq. (Q) and the width parameters o", for i = o,d as obtained by the 
fit ([17|) we determine finite volume estimators of the specific heats Co{L) 
and Cd{L). The results are collected in Table 3 for q=10 and in Table 4 
for q=20. Figure displays a plot of Co{L) for q=10 and q=20, together 
with asymptotic values predicted by a large q-expansion [0. In case of 
q=20 our data for Co{L) converge towards a value consistent with the large 
q expansion results and the FSS study |^ that gives Co = 5.2(2). In addition, 
the data for Cd{L) — Co{L) approach, as L grows, the exactly known infinite 
volume limit. For the q=10 model the data overshoot the value Co = 12.7(3), 
obtained from the FSS of the extrema of the specific heat |1T3. The value 



indicated by the large q-expansions is Co = 18.06(4). A value Co ~ 18 was 
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obtained in []T3[ from the FSS of the specific heat at Pf Much larger lattices 
would be required in order to see Co{L) approach this limit (The correlation 
length 1^ at (3t is ^disorder = 10.56 for g = 10 and 2.70 for q = 20). 

In figure ^ to |^ we display the locations of the maxima E^""^'" and El^"'^''^ 
of the energy probability distribution function at f5t as function of the variable 
1/L'^, compared with the prediction E^°'^'^ = Ei + const./ L'^ with slopes [|12 



100 (g = 10) and 15 (g = 20). The data are collected also in Tables 3 

and 4. For our largest lattices the predicted linear dependence on 

is born out. The deviations observed can be reproduced within the 



formalism of [Q. This rules out the 1/L behavior predicted in |]T5|. For 
the larger sized lattices the ratio Rq become consistent with the theoretical 



prediction q^^Cd/ Co- These tables contain also our estimates for the partition 
functions ratios Z°j^{(3t) / Zj^i^Pt) ■ They are defined by summing the probability 
distribution function at j3t over corresponding energy ranges E < i?™", resp. 
E > E^^". These ratios are consistent with the value q as predicted by 
theory. Unfortunately, due to small numbers of tunneling events, they have 
sizable errors on our largest lattices. 

In figure H we display a selection of our data for the probability distribu- 
tion function Pi(/3i,£') for the q=20 Potts model. An analogous figure for 
the q=10 Potts model can be found in 0. We observe the unfolding of a flat 
region of the probability distribution function on our largest lattices L > 38, 
therefore giving direct numerical evidence for the dominance of two planar 
interface configurations. On the 50^, 60^ and 70^ lattices we have tried to 
estimate the limits of the energy range where this occurs. For this purpose 
we have numerically estimated the curvature of the probability distribution 
function. We assume that on finite lattices the locations of largest curvature 
determine the limits. An example of this procedure is displayed in figure ^ 
The obtained values, as collected in Table 5, are rather close to the values 
E^ of the spherical bubble model. 

Tables 3 and 4 also contain our results for the finite volume estimators of 
the interface free energy F{L). They have been plotted in figures |10] and [Tl| 
as a function of the variable 1/L. The horizontal line in these figures always 
denotes the known analytical value [J^. It can be seen that in case of the 
q=10 model the finite size behavior of F{L) is very well described, in the 
whole L range, by a correction proportional to 1/L {i.e. p = 1) as predicted 
by the capillary wave approximation. However in the q=20 Potts model 
we find a more complicated behavior. We cannot draw a conclusion on the 
actual form of the functions F{L), i.e. it is undecided whether logarithmic 
corrections are present or not. 

We have performed several fits with the form eq. (p!6D . We have succes- 
sively freed higher order corrections in 1/L and also varied the interval of 
lattices sizes. The results of this fits together with their /d.o.f. values are 
collected in Tables 6 and 7. All fits of Table 6 and most fits of Table 7 are 
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within one standard deviation consistent with the exact Fo,d. result of table 
1 (In the g = 10 case, this consistency is obtained although we do not ob- 
serve the expected flat part in PL^Pt, E), in the g = 20 case this consistency 
requires either to omit small lattice data, or to include corrections). 
As our final result we quote 

Fo.d = 0.0950(5) for g = 10 (18) 

and 

Fo.d = 0.3714(13) for q = 20. (19) 



5 Conclusion 

We have presented a careful analysis of properties of the probability distribu- 
tion function Pl(/3, E) in the 2-d q=10 and q=20 Potts model at the infinite 
volume transition point. 

• On our largest lattices the ratio of the heights of the two maxima 
of PL{Pt,E) follows the prediction of the Gaussian model. The ratio of 
the ordered to the disordered partition function is equal to q. The finite 
volume corrections to the locations of the maxima of the energy follow a 
l/L*^ behavior with a slope as predicted by the large q expansion. 

• The Gaussian width of the peaks of PidSt, E) approaches for large 
volumes the prediction of the Gaussian model with specific heat parameters 
Co and Cd, which in the case of the q=20 Potts model are very well consistent 
with the large q expansion and a determination using finite size scaling. In 
the g = 10 case lattices much larger then 100^ would be required in order to 
see the asymptotic behavior. 

• The finite volume interface free energy estimators allow in both cases 
a natural extrapolation of the data to the recently calculated exact infinite 
volume values. Our results are consistent with those values. 

• Around iJ™" we observe the unfolding of a fiat region in PiiPt, E). It 
is an interesting recent observation that this effect can be enhanced by using 
lattices which are elongated in one direction . 
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q 


Fo.d. 


MC estimate 


7 


0.020792 


0.0241 (10) fig 


10 


0.094701 


0.0978 (08) § 


20 


0.370988 





Table 1: The exact interface free energy for the 2-d q=7, 10 and 20 Potts 
models, together with previous MC estimates. 



L 


Msweeps 


50 


26 


60 


22 


70 


10 



Table 2: Additional statistics on 50^, 60^, and 70^ lattices in the q=20 Potts 
models in units of million sweeps. 



L 


Co 

Cd 


Cd — Co 


j^max,d 
j^max,o 


-Rio 

Zo/Zd 


F{L) 


12 


4.81(06) 
5.77(09) 


0.97(12) 


-.8813(08) 
-1.7666(05) 


10.70(10) 
11.00(09) 


0.1379(04) 


16 


5.33(04) 
6.45(06) 


1.09(08) 


-.9012(04) 
-1.7399(04) 


10.25(15) 
9.76(13) 


0.1273(04) 


24 


7.06(05) 
7.84(06) 


0.77(09) 


-.9243(06) 
-1.7112(05) 


9.98(28) 
9.46(25) 


0.1168(05) 


34 


8.64(04) 
9.43(05) 


0.82(06) 


-.9391(02) 
-1.6954(02) 


10.39(57) 
9.88(51) 


0.1098(07) 


50 


10.66(09) 
11.14(11) 


0.37(16) 


-.9504(03) 
-1.6835(03) 


10.56(74) 
10.33(69) 


0.1056(06) 


70 


12.60(14) 
13.08(23) 


0.48(28) 


-.9568(05) 
-1.6760(04) 


10.2 (12) 
9.7 (11) 


0.1022(09) 


100 


13.99(20) 
14.71(35) 


0.54(45) 


-.9612(05) 
-1.6712(02) 


9.2 (17) 
9.1 (16) 


0.0997(11) 


oo 


18.06(4) 
18.51(4) 


0.44763 


-.96820 
-1.66425 


10.12 
10.0 


0.094701 



Table 3: Finite volume results for the q=10 Potts model at (3t, together with 
theoretical expectations in the infinite volume limit. 
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L 


Co 

Cd 


Cd — Co 


^max,d 

j^max,o 


R20 

Zo/Zd 


m 


16 


3.473(10) 
4.215(08) 


0.742(14) 


-.60079(09) 
-1.85270(09) 


22.40(16) 
20.08(13) 


0.3772(02) 


18 


3.486(09) 
4.403(12) 


0.912(16) 


-.60482(13) 
-1.84726(10) 


22.71(28) 
20.41(25) 


0.3762(03) 


20 


3.768(32) 

4.556(10) 


0.789(36) 


-.60796(09) 

-1.84307(10) 


22.19(26) 

20.16(23) 


0.3762(03) 


24 


3.895(08) 
4.808(08) 


0.909(12) 


-.61236(07) 
-1.83765(06) 


22.62(24) 
20.61(22) 


0.3759(02) 


30 


4.196(12) 
5.109(12) 


0.908(16) 


-.61653(09) 

-1.83265(07) 


21.44(39) 
19.66(36) 


0.3751(03) 


32 


4.270(13) 
5.213(14) 


0.942(17) 


-.61761(09) 
-1.83149(07) 


21.66(65) 
19.85(59) 


0.3758(03) 


34 


4.370(26) 
5.204(36) 


0.831(35) 


-.61856(14) 
-1.83048(11) 


22.87(46) 
20.98(36) 


0.3766(10) 


38 


4.495(13) 
5.356(18) 


0.863(24) 


-.61980(08) 
-1.82877(06) 


20.84(68) 
19.24(63) 


0.3747(05) 


50 


4.739(20) 
5.603(26) 


0.864(36) 


-.62236(14) 
-1.82572(09) 


23.9 (24) 
22.1 (22) 


0.3747(09) 


60 


4.925(19) 
5.693(30) 


0.768(30) 


-.62334(14) 
-1.82410(10) 


26.0 (41) 
24.3 (38) 


0.3738(09) 


70 


5.004(28) 

5.825(31) 


0.821(42) 


-.62430(13) 

-1.82334(10) 


27.0 (45) 

25.1 (42) 


0.3732(13) 


oo 


5.362(3) 
6.133(3) 


0.77139 


-.62653 
-1.82068 


21.5 
20.0 


0.370988 



Tabic 4: Finite volume results for the q=20 Potts model at (3t> together with 
theoretical expectations in the infinite volume limit. 



L 


Et 


EI 


50 


-1.00 


-1.35 


60 


-1.00 


-1.32 


70 


-1.01 


-1.36 


bubble 


-1.01 


-1.44 



Tabic 5: Prediction of the spherical bubble model (last line) together with 
estimated locations on 50^, 60^ and 70^ lattices in the q—20 Potts model. 



12 



nt 


r 

-L-'min 


^o.d. 


^ ,2 / J „ f 


a2 — as — — u[ ) 




u.uyoui^uo ) 


.ZLo 




16 


0.0949(06) 


.220 




24 


0.0946(08) 


.206 




34 


0.0949(13) 


.257 




50 


0.0938(20) 


.003 


aa = 6 = 


12 


0.0946(09) 


.176 




16 


0.0943(13) 


.218 




24 


0.0947(23) 


.307 




34 


0.0924(39) 


.041 


6 = 


12 


0.0942(19) 


.221 




16 


0.0945(33) 


.325 




21 


0.0905(61) 


.103 


02 = 03 = 


12 


0.0941(15) 


.167 




16 


0.0939(22) 


.217 




24 


0.0945(40) 


.308 




34 


0.0903(69) 


.062 


Theory 




0.094701 





Table 6: Results of fits of the form eq.(16) to the finite volume estimators 
F{L) of Table 3 for the q—10 Potts model. Column one specifies which pa- 
rameters of the fit where fixed to 0, with b — —p + d/2. The fit intervals are 
from L„iin of column two to Lmax = 100. The fit marked with (*) denotes the 
fit which results into our infinite volume estimate eq.(18). It has been plotted 
in figure 10. 
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fit 




Fo.d. 


x'/d.o.f 


consistent 
with theory 


02 = 03 = ^' = 


16 


0.3736(04) 


1.187 


no 




18 


0.3739(05) 


1.184 


no 




20 


0.3736(06) 


1.179 


no 




24 


0.3733(07) 


1.306 


no 




30 


0.3726(12) 


1.470 


no 


(*) 


32 


0.3714(13) 


.590 


yes 




34 


0.3713(17) 


.782 


yes 




38 


0.3724(19) 


.265 


yes 




50 


0.3694(45) 


.001 


yes 


as = b = 


16 


0.3733(12) 


1.330 


no 




18 


0.3717(15) 


.985 


yes 




20 


0.3716(17) 


1.149 


yes 




24 


0.3709(23) 


1.329 


yes 




30 


0.3663(46) 


1.335 


no 




32 


0.3726(59) 


.773 


yes 




34 


0.3766(94) 


1.009 


yes 




38 


0.362 (13) 


.028 


yes 


6 = 0(**) 


16 


0.3684(29) 


1.033 


yes 




18 


0.3701(38) 


1.116 


yes 




20 


0.3687(46) 


1.287 


yes 




24 


0.363 (10) 


1.528 


yes 




30 


0.390 (21) 


1.368 


yes 




32 


0.358 (31) 


1.057 


yes 


02 = 03 = 


16 


0.3727(20) 


1.312 


yes 




18 


0.3698(25) 


.962 


yes 




20 


0.3699(30) 


1.121 


yes 




24 


0.3688(40) 


1.312 


yes 




30 


0.3608(88) 


1.383 


no 




32 


0.373 (11) 


.777 


yes 




34 


0.379 (18) 


1.064 


yes 




38 


0.354 (25) 


.037 


yes 


Theory 




0.370988 







Table 7: Results of fits of the form eq.(16) to the finite volume estimators 
F{L) of Table 4 for the q=20 Potts model. Column one specifies which pa- 
rameters of the fit where fixed to 0, with b — —p + d/2. The fit intervals are 
from Ljnin of column two to L^ax = 70. The last column specifies whether the 
exact Fo,d. result is within one standard deviation of the fitted value. The fit 
marked with (*) denotes the fit which results into our infinite volume estimate 
eq.(19), the one marked with (**) has been plotted in figure 11. 
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E 

FIG . 1 

Figure 1: Plot of In PiiPt-i E) on a 50^ lattice for the q=20 Potts model in 
the vicinity of the disordered state. The crosses denote (selected) data with 
their error bars. The circle denotes the maximum value P'^'^^'^ with error 
bar. The two curves display a cubic polynomial fit to In PiiPt, E) and the 
Gaussian contribution to this fit. 



15 




-1,2 -1,18 -1,16 -1,14 -1,1 



E 

FIG , 2 

Figure 2: Plot of In PL{/3t, E) on a 50^ lattice in the q=20 Potts model 
in vicinity of p™*". The crosses denote selected data. The wiggling curves 
indicate the error interval of the data, while the circle with errror bars denotes 
our estimate of P™" and ii;™*". 
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A q=10 Potts model 
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q=20 Potts model 
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FIG, 3 



Figure 3: Estimates of Co for the q=10 (circles) and q=20 Potts (triangles) 
models, as extracted from the width of the ordered peak of Phil^t, E). The 
full symbols denote asymptotic results from the large q-expansion. 
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008 



FIG , 4 



Figure 4: Locations of the maxima of disorderd states at A in the 

q=10 Potts model as function of the inverse volume. The full symbol and 
dashed line denote the infinite volume prediction. The full line has a slope 
g'" /{2g") as predicted by the large q expansion [12]. 
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FIG, 5 



Figure 5: Locations of the maxima of ordered states E^"'^'" at (3t in the 
q=10 Potts model as function of the inverse volume. The symbols and hues 
have the same meaning as in Fig. 4. 
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Figure 6: Locations of the maxima of disordered states E^^"'^''^ at /3j in the 
q=20 Potts model as function of the inverse volume. The symbols and hues 
have the same meaning as in Fig. 4. 
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Figure 7: Locations of the maxima of ordered states E^"'^'° at (3t in the 
q=20 Potts model as function of the inverse volume. The symbols and hues 
have the same meaning as in Fig. 4. 
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FIG. 8 

Figure 8: Plot of In PL{(3t,E) for selected lattice sizes in the q=20 Potts 
model. The unfolding of a fiat region with increasing lattice size L > 38 is 
observed. 
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FIG , 9 



Figure 9: Plot of In PL{f3t,E) on a 70^ lattice. The statistical errors cannot 
be resolved. In the upper half of the plot we display a numerical estimate of 
the curvature of In E). The transition from the two planar interface 

to the single bubble region is indicated by two peaks. Their position is 
compared with the prediction for spherical bubbles . 
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infinite volume theory 
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1 



Figure 10: Plot of the finite volume estimator F{L) for the q=10 Potts model 
as function of 1/L. The horizontal hne corresponds to the theoretical value, 
and a fit of the form F(L) — F^d. + ai/L is displayed. 



24 




Figure 11: Plot of the finite volume estimator F{L) for the q=20 Potts model 
as function of 1/L. The horizontal hne corresponds to the theoretical value, 
and a fit of the form F(L) — F^d. + ai/L + 02/-^^ + a^/ is displayed. 
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